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Abstract 


We have examined cutoffs and pile-ups due to various processes in the spectra of 
particles produced by shock acceleration, and found that, even in the absence of energy 
losses, the shape of the spectrum of accelerated particles at energies well below the nominal 
maximum energy depends strongly on the energy dependence of the diffusion coefficient. 
This has implications in many areas, for example, in htting the observed cosmic ray 
spectrum with models based on power-law source spectra and rigidity dependent diffusive 
escape from the galaxy. With continuous energy losses, prominent pile-ups may arise, 
and these should be included when modelling synchrotron X-ray and inverse Compton 
gamma-ray spectra from a shock-accelerated electron population. 

We have developed a Monte Carlo/numerical technique to model the shape of the 
spectrum for the case of non-continuous energy losses such as inverse Compton scatter¬ 
ing in the Klein-Nishina regime. We hnd that the shapes of the resulting cut-offs differ 
substantially from those arising from continuous processes, and we suggest that such dif¬ 
ferences could be observable through their effect on the spectrum of radiation emitted by 
a population of recently accelerated electrons as, for example, may exist in young super¬ 
nova remnants. 

PACS codes/keywords: 96.40 Cosmic rays; 96.40.De Composition, energy spectra, 
and interactions; 98.70.Sa Cosmic rays (including sources, origin, acceleration, and inter¬ 
actions) 
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1 Introduction 


Cosmic rays with energies up to 100 TeV are thought to arise predominantly through shock 
acceleration by supernova remnants (SNR) in our Galaxy [|^. A fraction of the cosmic rays 
accelerated should interact within the supernova remnant and produce gamma-rays P,|^. 
The EGRET instrument on the Compton Gamma Ray Observatory has observed gamma 
ray signals above 100 MeV associated with at least two supernova remnants - IC 443 
and 7 Cygni |^. (It is however possible that the gamma ray emission from IC 443 is 
associated with a pulsar within the remnant rather than the remnant itself [Q). Further 
evidence for acceleration in SNR comes from the recent ASCA observation of two spatially 
resolved regions with non-thermal X-ray emission from SN 1006 These non-thermal 
X-rays are interpreted as synchrotron emission by electrons accelerated in the remnant 
up to energies as high as 100 TeV [0, |[. Donea and Biermann ^ however suggest it 
may be bremsstrahlung from much lower energy electrons. The spectrum of synchrotron 
and inverse Compton emission will depend on the the electron spectrum. As we shall 
demonstrate in this paper the effect of a cut-off in the spectrum due to escape losses, 
energy losses or interactions can considerably modify the shape of the electron spectrum 
over several decades in energy below the nominal cut-off energy. 

In this paper, we adopt a simple leaky-box model of shock acceleration based on the 
physical picture of acceleration at plane non-relativistic shocks. This simple model will 
then be used to explore the effect of cut-offs due to the hnite size of the accelerator and 
energy-loss and interaction processes. 


2 Leaky-box picture of shock acceleration 


The leaky-box picture of shock acceleration |]^ is based on a simple heuristic treatment 
of Fermi acceleration. See for a simple treatment based on that in Gaisser’s book [^. 
More detailed and rigorous treatments are given in several review articles |]I4|-[T7[ . 

The rate of gain of energy of relativistic particles at non-relativistic shocks is given by 


dE 

dt 


{AE) 


( 1 ) 


where Cycle is the time for one complete cycle, i.e. from crossing the shock from up¬ 
stream to downstream, diffusing back towards the shock and crossing from downstream 
to upstream, and hnally returning to the shock, and (AE) is the average energy gain per 
cycle. 

By considering the collisionless scattering of relativistic particles by magnetic scatter¬ 
ing centres in the upstream and downstream plasmas, one hnds 


{AE) ^ 4{R-l) Vs 
E 3R c 


2 






















where Vs is the shock velocity, R = ui/u 2 is the compression ratio, ui = Vs and U 2 are the 
upstream and downstream velocities in the rest frame of the shock. For a strong shock in 
a monatomic gas (e.g. fully ionized plasma) R = 4. The acceleration rate at energy E is 
dehned by 


- 

E dt 


(3) 


The rate of loss of accelerated particles downstream is the probability of escape per 
shock crossing divided by the cycle time 


r 


esc 


Prob. (escape) 


t 


cycle 


(4) 


By considering the rate at which cosmic rays cross the shock from upstream to down¬ 
stream, and rate of convection of cosmic rays downstream one obtains the probability of 
escape downstream per shock crossing 


Prob. (escape) ~ 


R c 


(5) 


We see immediately that the ratio of the escape rate to the acceleration rate depends 
on the compression ratio 

^esc 3 , , 

-~ -B - 7 ( 6 ) 

^ a.cc 1 


and for a strong shock {R = 4) the two rates are equal. As we shall see below, a conse¬ 
quence of this is that the asymptotic spectrum of particles accelerated by a strong shock 
is the well-known E~‘^ power-law. 

The cycle time depends on the diffusion coefficients upstream, ki, and downstream, 
/c 2 , of the shock, and on the shock velocity and compression ratio and is given by 



The acceleration rate is then given by 


{R - l)ui f h 
3R U2 


-1 


(7) 

( 8 ) 


Assuming the the diffusion coefficients upstream and downstream have the same power- 
law dependence on energy, e.g., 

ki (X k2 (X E^, (9) 

then the acceleration rate also has a power-law dependence 

race oc E~\ (10) 


Note that more correctly the diffusion coefficient will be a function of magnetic rigidity, 
p, which, for ultra-relativistic particles considered in this paper, is approximately equal 
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to E/Ze where Ze is the charge. However, in this paper since we are mainly concerned 
with singly charged particles we shall work in terms of E rather than rigidity. 

The leaky-box acceleration model may then be considered as follows. A particle of 
energy Eq is injected into the leaky box. While inside the box, the particle’s energy 
changes at a rate dE/dt = Era,cc{E) and that in any short time interval At the particle 
has a probability of escaping from the box given by Atresc{E). The energy spectrum 
of particles escaping from the box then approximates the spectrum of shock accelerated 
particles. 

2.1 Energy spectrum for acceleration rate different to escape 
rate 

Let us consider hrst the case of no energy losses, or losses due to any other process. Nq 
particles of energy Eq are injected at time t = 0, and we assume the following acceleration 
and escape rates: 

race = aE~\ (11) 

Tesc = cE~^. ( 12 ) 

The energy at time t is then obtained simply by integrating 

dE/dt = aE^^-^\ (13) 

giving 

E{t) = {E/,+6aty/^. (14) 

The number of particles remaining inside the accelerator at time t after injection is ob¬ 
tained by solving 

dN/dt = -N{t)c[E{t)]-\ (15) 

Using Eq. [1^ and integrating, one has 

rN{t) pt 

/ N-^dN = -c / {El + 5at)-^dt, (16) 

J Nq J 0 

giving 

N{t) = No[E{t)/Eo]-^/\ (17) 

Since Nq — N{t) particles have escaped from the accelerator before time t, and therefore 
have energies between Eq and E(f), the differential energy spectrum of particles which 
have escaped from the accelerator is simply given by 

dN/dE = Nq{T-1){Eq)-\E/Eq)-^, {E>Eq) (18) 

where T = (1 -f c/a) is the differential spectral index. We note that for resc{E) = ra,cc{E) 
one obtains the standard result for acceleration at strong shocks T = 2. 
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3 Cut-off due to finite acceleration volume, etc. 

Even in the absence of energy losses, acceleration usually ceases at some energy due to the 
hnite size of the acceleration volume (e.g. when the gyroradius becomes comparable to 
the characteristic size of the shock), or as a result of some other process. We approximate 
the effect of this by introducing a constant term to the expression for the escape rate: 

resc = cE-^ + (19) 


where is dehned by the above equation and will be close to the energy at which 

the spectrum steepens due to the constant escape term. We shall refer to as the 
“maximum energy” even though some particles will be accelerated to energies above this. 

Following the same procedure as for the case of a purely power-law dependence of 
the escape rate, we obtain the differential energy spectrum of particles escaping from the 
accelerator. 


dN 

Je 


= iVo(r-l)(Eo)-^(E/Eo) 
x[l -h (E/Emax)'^] exp 


\-r 


r-1 


E 


E„ 


En 


E„ 


[E > Eq) (20) 


for 6 > 0. For (5 = 0 we note that from Eq. = (2c) at all energies. Thus the 

situation is equivalent to the case of F^ax —^ oo and the spectrum is given by Eq. |T8 
provided we replace c with 2c, i.e. we must replace T with T' = ( 2 r — 1). 

We show in Fig |I| the differential energy spectrum for T ranging from 1.5 to 2.5, for 
(a) 6 = 1/3, (b) 6 = 2/3, (c) 5 = 1. In Fig ^ we compare the spectra for T = 2 and 6 
ranging from 1/3 to 1, and note that the energy dependence of the diffusion coefficient 
has a profound influence on the shape of the cut-off. 

Spectra such as those presented in Fig. [| may be used to model the source spectra 
of high energy cosmic ray nuclei of various species if one replaces energy with magnetic 
rigidity, p = pc/Ze ~ F/Ze, where Z is the atomic number: 


dN 

dp 


= iv„(r-i)(p„)-'(p/p(,) 


[1 + (P/Pmax)'^] exp 


\-r 


r-1 


p 


Po 




Pn 


ip>Po) (21) 


4 Cut-off due to E‘^ energy losses 

E"^ energy losses of electrons result from synchrotron radiation and inverse Compton 
scattering of low-energy photons in the Thomson regime. This may be treated as a 
continuous energy loss process. At higher energies, when the electron rest frame energy 
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of the emitted photon is no longer much less than the electron rest mass energy (Klein- 
Nishina regime of Compton scattering, or quantum synchrotron regime), the energy losses 
are no longer continuous and this will be discussed in Section 5. 

The introduction of continuous energy losses into the problem is, in principle, straight¬ 
forward and accomplished simply by modifying Eq. p!3| , 

dE/dt = -bE^. (22) 


For this case, however, the problem is easier to solve using the Green’s function approach 
adopted by Stecker [^] when considering the ambient spectrum of cosmic ray electrons 
and galactic gamma-rays. Using the appropriate Green’s function one can obtain the 
steady-state spectrum of particles inside the leaky-box accelerator, and multiplying this 
by the escape rate one obtains the spectrum of particles leaving the accelerator. 


where 


dN 

He 


cE ^ + E^ 


-5 

max 


',E^ 


-5 


bE^ 


exp[-/(E)] 


HE) ^ [ 

Je 


® c E-^ + E, 


^-S 


Eo aE^-^-bE^ 


dE 


(23) 


(24) 


is given in Appendix A. The result depends on the parameters 6, T, Eq, Ecut, 

As a result of the energy loss by particles near the maximum energy a pile-up in the 
spectrum may be produced just below E^nt- The size of the pile-up will be determined by 
the relative importance of race and r^sc at energies just below Ecut, and so should depend 
only on T and 5 provided Eq -C Ecut -C E mav In this case, the shape of the spectrum is 
given by 


dN 

He 


]v„(r-i)(B„)-‘(|-) 


-r 


^ (r-2-5)/(i+5) 


E 


cut 


(25) 


We compare in Fig. the spectra for this case for 5 = 1/3, 2/3 and 1, and for T = 1.5, 
2, and 2.5, and note that the pile-ups are higher for flatter spectra, and that for steep 
spectra the pile-up may be absent or the spectrum may steepen before the cut-off if 5 is 
small. The effect of synchrotron losses on multiple diffusive shock acceleration have been 
considered by Melrose and Grouch [^, and similar effects were found. 

For the case where Ecut is comparable to or not much less than Emax we show in Fig ^ 
the differential energy spectrum for T = 2 for (a) 5 = 1/3, (b) 5 = 2/3, (c) <5=1, and 
in each case we show results for Ecut/Emax ranging from 10“® to 10^. Figs. H(a)— (c) 
are identical in the sense that they are all for T = 2 and show synchrotron cutoffs at 1, 
1/10, ... 10“® of E m a x. The only difference is in the energy dependence of the diffusion 
coefficient. In Fig. H(c), for 5=1, the shape of the pile-up is independent of Ecut for 
Ecut < O.lEjuax- In contrast, in Fig. |Ka), for 5 = 1/3, the size of the pile-up decreases 
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with increasing h^cutZ-f^max, and so does not depend simply on T and S. This is because 
the constant escape term gives rise to a steepening in the spectrum, and thus from Fig. ^ 
we expect a smaller synchrotron pile-up, or even a steepening. For example, note that 
in Fig. §(a), for 5 = 1/3 the the spectrum in the absence of synchrotron losses steepens 
at signihcantly lower energies (due to the constant leakage term) than for the case of 
5 = 1, and a.t E = 10“^i?max has a differential index of ~ 2.2. We plot in Fig. |^(d) results 
generated for F = 2.2 and 5=1/3 and note that the size of the pile-up for Ecnt ^ -^max 
is similar to the case for F = 2, 5 = 1/3, and E = 10“^i?max- 


5 Cut-off due to inverse Compton losses on the 2.7 K 
background 

In this section, we discuss the acceleration of electrons in the presence of a blackbody 
radiation held. Such a scenario is important in many areas of high energy astrophysics. For 
example, one may have acceleration of electrons in active galactic nuclei being responsible 
for radio, optical. X-ray, and gamma-ray emission, and this acceleration takes place in 
the presence of an intense radiation held, e.g. from the accretion disk. 

We choose as our example, the case where the cosmic microwave background radiation 
is the dominant radiation held. Such a situation would occur during 1st order Fermi 
acceleration of cosmic ray electrons in our Galaxy at a shock in the interstellar medium 
due to a young supernova remnant. 


5.1 Times scales for losses and acceleration 


The energy loss rate is given by 

f = -rrcU,^ ( 26 ) 

where ar is the Thomson cross section, 7 is the electron’s Lorentz factor, and U = /'^tt 

for synchrotron losses and U is the energy density of the radiation held for inverse Compton 
losses in the Thomson regime (Grad = for blackbody radiation). Thus, at low energies, 
the time scale for energy loss of an electron by synchrotron losses is given by 


= 1.3 X 10 


10 


—dE/dt 

while that for IC scattering is given by 


IGeVj llpG 


-2 


years. 


E 


-dE/dt 


= 1.3 X 10^ 


-^y 

1 GeVy 


T 


2.735 


-4 


years. 


(27) 


(28) 


The energy loss time scale for synchrotron losses is shown in Figs. |^(a)— (c) for B ranging 
from 10“^^ G to 10“^ G by the dotted lines. 
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For electrons with energies above ~ 3x10^ GeV energy losses during inverse Compton 
scattering on the microwave background photons are no longer continuous. The Klein- 
Nishina cross section must be used, and in each interaction the electron loses a substantial 
fraction of its energy. Monte Carlo simulations using the differential Klein-Nishina cross 
section have been performed and the mean energy loss per interaction (AE) has been 
calculated as a function of the electron energy. One can thus obtain an effective energy- 
loss time scale, Et[^t/(AE), where tint is the mean interaction time. See [p5|-^ for 


discussion of Monte Carlo calculations of inverse Compton scattering. The energy loss 
time scale and the interaction time for IC losses are shown in Figs. |^(a)— (c) by the thick 
solid curves, and dashed curves labelled TC’ (the thin solid curves, and dashed curves 
labelled ‘TPP’ give the energy loss time scale and the interaction time for triplet pair 
production; see e.g. [p^). 

During acceleration, the cut-off energy (and the loss process responsible) will depend 
on the strength of the magnetic held and the acceleration rate. The acceleration rate is in¬ 
versely proportional to the diffusion coefficient (see Eq. |^), and so the highest acceleration 
rate is achieved with the lowest diffusion coefficient. 

The diffusion coefficients required ki and /c 2 are the coefficients for diffusion parallel 
to the shock normal. The diffusion coefficient along the magnetic held direction is some 
factor Tj times the minimum dihusion coefficient, known as the Bohm dihusion coefficient. 


*•'11 = ’q-Rpc 


(29) 


where Rg is the gyroradius, and rj > 1. 

Parallel shocks are dehned such that the shock normal is parallel to the magnetic held 
{B\\ui). In this case, making the approximation that ki = k2 = fc||, i? = 4, and Bi = B2 
one obtains 


3 / 9 ? 

20ri^^ 


(30) 


where jSi = U\/c and ojg 
an acceleration rate 


c/Rg. For a shock speed of Ui 

do « 1-5 X lO-L..,, 


0 .1c and rj = 10 one obtains 


(31) 


For the oblique case, the angle between the magnetic held direction and the shock 
normal is diherent in the upstream and downstream regions, and the direction of the 
plasma how also changes at the shock. The dihusion coefficient in the direction at angle 
6 to the magnetic held direction is given by 


k = kn cos^ 9 + k 2 _ sin^ 9 


(32) 


where k± is the dihusion coefficient perpendicular to the magnetic held. Jokipii shows 
that 


k 


1 + rf 


(33) 














provided that rj is not too large (values in the range up to 10 appear appropriate). 

In the case of acceleration at perpendicular shocks, Jokipii [^] has shown that accel¬ 
eration can be much faster than for the parallel case. Using k± and B 2 ~ 4i?i and one 
obtains 

... 
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For a shock speed of Ui = 0.1c and rj = 10 one obtains an acceleration rate of 


^acc ~ 


O.OAuJn 


(34) 


(35) 


Generally, for stochastic particle acceleration by electric helds induced by motion of 
magnetic helds the acceleration rate is given by 


r-acc = i{E)ujg ^ O.Ol^(U) 


-1 


B 


1 /iG 


(36) 


where ^{E) < 1 and depends on the acceleration mechanism; constant values of ^ ranging 
from 1.5 X 10“^ to 0.04 might be achieved by hrst order Fermi acceleration with Bohm 
diffusion [k cx U, i.e. 5 = 1). 

We show in Fig. ^(a), along with the energy loss time scales for IG and synchrotron 
losses, the acceleration time scale for 5=1 (Bohm diffusion) for ranges of ^B which are 
probably appropriate for supernova shocks. Gomparing the energy loss time scale for IG 
and synchrotron loss we see that if i? > 3 /iG the cut off will be determined by synchrotron 
losses irrespective of the acceleration rate. Gonversely, if i? < 3 /iG the cut off will be 
determined by IG provided the acceleration rate is low such that the acceleration time 
curve intersects the IG loss curve in the Thomson regime. However, if the magnetic held 
is less than ~ 1 /iG and ^B is in the range 10“® to 10“® /xG, the cut-off in the spectrum 
of accelerated electrons will be determined by IG scattering in the Klein-Nishina regime. 
If the magnetic held were less than 0.1/xG, the range of ^B would be extended up to 
10“® /xG, but for ^B > 1.4 x 10“^/xG the cut oh will be determined by synchrotron losses. 
If the cut oh is be determined by IG in the K-N regime, the shape of the cut-oh would 
be diherent from that calculated for the case of E'^ energy losses, and one would need 
to take account properly of inverse Gompton collisions in the Klein-Nishina regime. The 
spectrum of synchrotron radiation and inverse Gompton gamma-rays would also depend 
on the details of the electron spectrum cut-oh and would therefore diher from the 
energy loss case. 

In Fig. |^(b) we show acceleration time scales that would apply if we used the dihu- 
sion coefficient inferred from a comparison of cosmic ray data with propagation model 
calculations 13 , k = 2.5 X 10^®(U/3 GeV)*^'® cm^ s and shock velocities ranging from 


0.01c to 0.1c. In this case, the cut-oh energy would probably be determined by IG in the 
Thomson regime or by synchrotron losses. 

In Fig. 11(c) we show acceleration time scales that would arise in the case of a Kol¬ 
mogorov spectrum of turbulence. A range of acceleration times (chain lines) are shown 
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for the purpose of illustration only. The diffusion coefficient can never be lower than the 
minimum (Bohm) diffusion coefficient, and so with increasing energy these curves would 
violate the Bohm limit. Thus the curves shown above ~ 10^ GeV are unrealistic for accel¬ 
eration at supernova remnant shocks in our Galaxy with typical galactic magnetic helds. 
The curves in Fig. ^c) should bend upwards and smoothly join an appropriate curve 
from Fig. ^(a). However, with an acceleration rate such as given by the middle curve, the 
cut-off energy determined by IG scattering would not be precisely determined, and could 
be anywhere between 3 x 10^ and 3 x 10^ GeV, and this would result in an interesting 
spectral shape. While this scenario does not apply at galactic supernova shocks, it may 
well apply in other astrophysical environments where shocks are present and different 
radiation helds exist. 


6 Monte Carlo Simulation 

A set of programs was constructed that simulates the particle acceleration in the leaky 
box acceleration model using the Monte Carlo method. In the simulation, a particle is 
injected with energy Eq and the acceleration and escape are simulated for a given F^max 
value and specihc energy dependence of the diffusion coefficient. After the injection with 
energy Eq the development of the electron energy spectrum is followed in time steps At. 
At every step the particles gain energy At{dE/dt). At was chosen to be small, 1/100 of 
the escape time at Eq. The spectrum of the escaping particles is calculated by weighting 
the particle energy after each step i by x P*“^, where P*“^ is the probability that the 
particle has remained in the accelerator region for i — 1 time steps. Pesc is calculated as a 
function of the particle energy for the chosen energy dependence of the diffusion coefficient 
with a constant Pesc(-Emax) added to simulate a smooth cutoff in the escaping particle 
spectrum. The maximum acceleration time was calculated so that particles could have 
reached energy of lOPmax in the absence of the constant term in the escape probability. 

The advantage of the Monte Garlo/numerical approach is that it is easily adapted 
to calculate the escaping electron spectrum for both continuous energy losses (e.g. syn¬ 
chrotron or IG in the Thomson regime) and non-continuous energy losses (e.g. IG in the 
KN regime). We next present results obtained by the numerical/Monte carlo technique 
produced in a series of runs made including synchrotron energy loss. The energy loss was 
calculated at every At and subtracted from the energy gain. The magnetic held value 
necessary for generating a cutoff at 1, 1/10, 1/100, 1/1000, and 1/10000 of Pmax was 
calculated keeping all ‘acceleration’ parameters as before. The program was tested for 
the same parameters used in the analytical/numerical technique described earlier, and 
the agreement was found to be excellent. We show, for example, in Fig. ^ the spectrum 
of accelerated particles in the presence of synchrotron losses for 5=1/3 obtained using 
the Monte Garlo method and compare this with the analytic results obtained previously. 
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6.1 Comparison of cutoffs due to synchrotron and IC energy 
loss 

Next we used the same procedure to simulate non-continuous inverse Compton energy 
loss. IC in the Klein-Nishina regime was incorporated as follows: for every time step, 
At, the IC interaction time was sampled from an exponential distribution with mean 
interaction time, tint as a function of the electron energy (given by the long-dashed curve 
in Fig. If the sampled interaction time was less than At, an IC collision was simulated 
by the Monte Carlo method and the electron energy loss was subtracted from the energy 
gain due to acceleration in the time step, i^max was set to 10® GeV to enable us to study 
the non-continuos energy loss in the Klein-Nishina regime. All acceleration parameters 
were adjusted to produce an IC energy loss time equal to the acceleration time at the 
same F'cut values as in the synchrotron case above. 

Figs. 0(a) and |0(b) show the electron spectrum in the presence of IC energy loss for 
the case of 5 = 1/3 and 5 = 1 respectively. Note that for Ecnt = 1000 GeV, where IG 
is in the Thompson regime and the energy loss is quasi-continuous, the picture is similar 
to that in the synchrotron case (Figs. 0(a) and (c)). For Ecnt = 10® GeV, where IG is 
in the Klein-Nishina regime, the shape of the cutoff is qualitatively different and much 
smoother. 

Fig. 0(b) is the version of Fig. 0(a) for the case of diffusion coefficient oc E. In 
the case of E^^^ energy dependence of the diffusion coefficient the pile-ups are much 
smaller and almost do not exist for Ecnt/Emnx greater than 10“^. Pile-ups are significantly 
more pronounced when the diffusion coefficient is proportional to E. The shape of these 
pile-ups, especially for IG losses in the Klein-Nishina regime, is very different from the 
synchrotron loss case (Figs. 0(a) and (c)). 

Fig. 0 compares the cutoffs due to synchrotron energy loss to those of IG energy loss 
on a linear scale. We show this comparison with a histogram because histograms are a 
better presentation of the Monte Garlo approach than smooth curves. 

The dotted histogram in Fig. || shows the cut-offs and pile-ups that are produced 
by IG energy loss. For Ecnt/Ecm.^ of 10“^ the shape of the pile-up in the spectrum of 
accelerated electrons is not very different from the synchrotron loss case. The reason 
is, of course, that the IG energy loss of 1000 GeV electrons in the Thompson regime 
is almost continuous, very similar to the synchrotron loss. For higher ratios, when the 
IG collisions move deeper and deeper in the Klein-Nishina regime and the energy loss 
becomes a significant fraction of the total electron energy, pile-ups become very wide, 
while for Ecnt = P^max there is no pile-up. It is replaced by a spectral cut-off which is 
very different from the cut-offs due to synchrotron loss or to escape from the shock. 

7 Conclusion 

We have derived the shape of the spectrum of particles produced by 1st order Fermi 
acceleration near the maximum or cut-off energy when the maximum energy is caused 
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by a constant escape rate term. We found that, even in the absence of energy losses, 
the shape of the spectrum of accelerated particles at energies well below the nominal 
maximum energy depends strongly on the energy dependence of the diffusion coefficient. 
This has implications in many areas, for example, in htting the shape of the observed 
cosmic ray spectrum in the region of the “knee” with model spectra based on power-law 
acceleration spectra and rigidity dependent diffusive escape from the galaxy. 

In the case where the cut-off is due to continuous energy losses a prominent pile-up 
may occur just below the cut-off energy dehned as the energy at which the total rate of 
change of energy is zero, and the spectrum cuts off sharply at this energy. 

We have developed a Monte Carlo/numerical technique to model cut-offs caused by 
non-continuous energy losses (i.e. interactions in which the particle being accelerated loses 
a substantial fraction of its energy). The Monte Carlo technique is tested by comparisons 
with the analytic calculations in cases where possible. Using numerical calculations we 
hnd that in the case of discrete energy loss the pile-up, if present, is less pronounced, and 
the cut-off is smoother with particle energies extending beyond the cut-off energy. This 
may well have observational consequences as it will affect the spectrum of synchrotron 
radiation and inverse Compton radiation. 

Further applications of this technique will compare the radiation spectra due to syn¬ 
chrotron radiation and to IC energy loss by the accelerated electrons and to those that 
escape from the shock in more realistic astrophysical environments, related to the environ¬ 
ments in different galactic supernova remnants. The technique developed can also be used 
to study the shape of the energy spectra of protons accelerated to ultra high energy in 
the presence of energy loss, such as pair production and photoproduction on the ambient 
photon fields. In addition, the Monte Carlo technique developed can be used for studies 
of the shape of the energy spectra of the accelerated particles in the presence of energy 
loss in more complicated acceleration models, including non linear shock models where 
the effective compression ratio is a function of the rigidity of the accelerated particles. 
We also intend to apply this method to study the shape of the proton spectra accelerated 
to ultra high energy in the presence of pair production and photoproduction losses in 
interactions on the ambient photon helds. 
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Appendix A 


Using the Green’s function (e.g. [^) for the “leaky-box model” of cosmic ray propagation, 
more properly called the homogeneous model, one can hnd the steady-state solution for 
the energy spectrum of particles inside the leaky-box 


f dE' \ 

dE R(E) \ Je„ t(E')R(E') J 


(37) 


where R{E) is the rate of energy gain and t{E) is the escape time. We can use this in 
the present case of the leaky box acceleration model with, R{E) = {aE^~^ + bE"^) and 
t{E) = {cE~^ + cE~l^)~^. We require the spectrum of accelerated particles, i.e. those 
escaping from the leaky-box, and this is simply dNi^/dE divided by the escape time t{E): 


dN 

Ie 


cE ^ + cEj^^ 

aE^-^ -bE^ 


iVoexp[-/(E)] 


where 


HE) 


[E cE ^ + cEj,^ 

Eo aE^-^-bE^ ’ 

^\eJ 


a6j \E, 


a(l -h 5) 


1 + 5 


, 1,1 + 


5 bE^+H 

1 + 5’ a ) 


In b 


E^+s 


Eq 


(38) 


(39) 


(40) 


and 2 E 1 is the hypergeometric function. 

Substituting the integration limits, and using F = (1 + c/a), and E^nt = {a/b)^/'d+^) 
we obtain 


dN 

Ie 




(41) 


where 


HE) = 


r -1 

5 

2 E 1 
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Er. 


1 + 5 


, 1,1 + 


2 EI 
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1 + 5 
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1 + 5’ \E, 
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^CUt 


l+<5- 


1 + 5’ \E, 


En 


^cut 


1+5- 
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In 


^-( 1 + 5 ) _ t^-{1+5) 
-^cut 


pi—(l+'5) 
-^cut 


^—(1+5) 


(42) 
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\og{E/E^J 


Figure 1: Differential energy spectrum for (a) 5 = 1/3, (b) 5 = 2/3, (c) 5 = 1. In each 
case we show results for F = 1.5 (highest curve), 1.6, 1.7, ..., 2.5 (lowest curve). 
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Figure 2: Differential energy spectrum for F = 2 and S = 1/3 (solid curve), 2/3 (dotted 
curve) and 1 (dashed curve). 
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Figure 3: Differential energy spectra for Eq Ecut ^ -Fmax for <5 = 1/3 (solid curves), 2/3 
(dotted curves) and 1 (dashed curves), and F = 1.5 (upper curves), 2.0 (middle curves) 
and 2.5 (lower curves). 
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Figure 4: Differential energy energy spectrum for F = 2 for (a) 5 = 1/3, (b) <5 = 2/3, (c) 
5 = 1, and (d) for F = 2.2 and 5 = 1/3. In each case we show results for i?cut/-Fmax = 10“® 
(leftmost curve), 10“^, ..., 10^ (rightmost curve). 
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log(time scale / years) log(time scale / years) 




log(£'/GeV) 

Figure 5: The energy loss time scale for synchrotron losses are shown by thin dotted 
lines for i? = 10 /xG (lowest line), 1 /iG, 0.1 /iG, ... 10“^^ G (highest line). For IG, the 
energy loss time scale (thick solid line labeled IG) and mean interaction time (dashed 
line labeled IG) are shown; the corresponding curves labelled TPP are for triplet pair 
production. Acceleration time scales are shown for: (a) 5 = 1 (Bohm diffusion) with 
= 1.4 X 10“® pG (leftmost solid line), 1.4 x 10“®, ... 1.4 x 10“^ pG (rightmost solid 
line); (b) 5 = 0.6 which has been suggested for the propagation of Galactic cosmic rays 
and Ml = 0.1c (lowest chain line), 0.03c, and 0.01c (highest chain line); (c) 5 = 1/3 which 
would arise in the case of a Kolmogorov spectrum of turbulence - a range of acceleration 
times (chain lines) are shown for the purpose of illustration. 
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log(V^max) 


Figure 6: Spectrum of accelerated electrons obtained using the Monte Carlo technique 
(histograms) compared with the analytic result from Fig. 4(a) for i?max = 10® GeV and 
cutoff due to synchrotron energy loss of 1, 1/10, 1/100, 1/1000, and 1/10000 i^max- The 
diffusion coefficient is proportional to 
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(a) 6 = 1/3 
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Figure 7: Spectrum of accelerated electrons for -Emax = 10' 
energy loss at 1, 1/10, 1/100, and 1/1000 -Emax- The diffusic 
to (a) (b) E. 
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GeV and cutoff due to IC 
1 coefficient is proportional 












\og{E/E^,:) 


Figure 8: Spectrum of accelerated electrons for -Emax = 10® GeV and cutoffs due to 
synchrotron (solid histogram) and IC energy loss (dotted histogram). Spectrum, divided 
by that for no IC loss, is shown on a linear scale to emphasize the cutoffs. The diffusion 
coefficient is proportional to {F/Fq represents the ratio of the spectrum calculated 

with energy losses to the spectrum calculated without energy losses). 


23 














